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Abstract 



The discovery of a heavy Higgs boson with mass up to mn = 1 TeV at the 
CERN LHC is possible in the H W^W" lujj decay mode. The weak 
boson scattering signal and backgrounds from tijj and from 1^+jets production 
are analyzed with parton level Monte Carlo programs which are built on full 
tree level amplitudes for all subprocesses. The use of double jet tagging and 
the reconstruction of the W invariant mass reduce the combined backgrounds to 
the same level as the Higgs signal. A central mini-jet veto, which distinguishes 
the different gluon radiation patterns of the hard processes, further improves the 
signal to background ratio to about 2.5:1, with a signal cross section of 1 fb. The 
jet energy asymmetry of the W jj decay will give a clear signature of the 
longitudinal polarization of the Ws in the final event sample. 
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I. INTRODUCTION 



In the effort to determine the dynamics of the spontaneous breaking of the electroweak gauge 
symmetry, the discovery of the Higgs boson would be of prime importance. Detecting the Higgs 
boson is one of the biggest challenges for the CERN LHC , both for a perturbative scenario 
for the symmetry breaking sector, with a Higgs boson mass below the Z-boson pair production 
threshold, and also if some strong interaction dynamics should be responsible for SU{2) x f/(l) 
breaking In both cases small signal rates, due to small usable decay branching fractions 

and/or small production rates, or large Standard Model (SM) backgrounds will have to be faced. 

In order to isolate a Higgs signal one will have to utilize all its characteristics. In turn this 
requires a simulation of the expected SM backgrounds with a high degree of detail, in a region of 
phase space where little or no experimental input exists at present. This problem is particularly 
acute for the search of a very heavy Higgs boson, with a mass above ~ 600 GeV. Here one will 
want to search for a Higgs resonance in the scattering of longitudinal weak bosons, or, more 
generally, one will look for some structure in the invariant mass distribution of the produced 
weak boson pairs in electroweak processes of the type qiq2 — * q^qiViVL @]. 

Numerous studies over the past several years have indicated that for the weak boson scattering 
signal to be identifiable, it is necessary to tag one or possibly two of the forward jets which arise 
from the scattered quarks A second characteristic of the weak boson scattering process is 

the lack of color exchange between the two incident quarks, which distinguishes it from typical 
background processes which proceed via the t-channel exchange of color octet gluons. These 
different color structures are expected to lead to a rapidity gap signature for the signal, either in 
terms of soft hadrons, at low luminosity ||10| , [TT| or in terms of mini-jets [|T^ . 

The small branching ratios of purely leptonic decays of the produced weak bosons can be 
overcome by studying the semi-leptonic modes, e.g. H — > W^W~ — > Ivjj. Here large back- 
grounds form ly+jets production put a premium on good W^-mass reconstruction of the two 
decay jets, in a situation where the large W energy leads to a small separation of the two jets. 
An advantage of this decay mode is the observability of the W jj decay angular distribution 
which may allow a measurement of the longitudinal W polarization of the signal. 

Most of these points have been considered before. The detectability of the H — > W^W~ — 
Ivji signal with jet tagging techniques, for example, has been discussed in the ATLAS and CMS 
technical design reports However, these studies have been based on parton shower Monte 

Carlo programs and it is not clear how well these programs describe the high pt jets associated 
with the decaying W . Also the color coherence effects which are at the basis of a rapidity gap 
trigger cannot be expected to be modeled correctly in these analyses. 

In this paper we perform a complementary study, based on full QCD matrix elements of all 
subprocesses contributing to the signal and to the various backgrounds. We consider the signal 
process 

qiq2 qsq^W^W q^q^ lyjj (1) 

(and crossing related ones) with a double forward jet tag on the two scattered quarks, ga and 
^4. For the dominant VT+jets QCD background we thus need the QCD matrix elements for all 
subprocesses leading to + 4 jet events |T^|T^. Similarly, the potentially large tt — > hW^hW~ 



background needs to be simulated with two additional partons in the final state [|T^], in order to 



account for the two tagging jets. When studying the consequences of different color structures 
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on soft gluon radiation patterns, the 0{as) QCD corrections for the signal must be known as 



well ||18|. While the parton level Monte Carlo programs for the individual subprocesses have 
been available in the literature, we here perform a first study of the H W^W~ Ivjj mode 
with full QCD matrix elements for signal and background subprocesses. 

In Section || we present these calculational tools in some detail. For the discussion of gluon 
radiation patterns we employ the truncated shower approximation which is briefly described 
at the end of that section. The isolation of the H — > W^W~ Ivjj signal, with double 
forward jet tagging, but without considering the ly-mass reconstruction from the W — >■ jj 
decay is considered in Section |T|. Here the hadronic system arising from the W-dec&y will be 
considered as a single jet. The properties of this W-dec&y jet, its internal dijet structure, and the 
measurement of the W^-mass is the subject of Section |^ Here we also consider the measurement 
of the W polarization via the energy asymmetry of the two decay jets [0. Parameterizing the 
results of the VT-mass analysis in terms of a reconstruction efficiency, we return to the simpler 
analysis, without simulating the W jj dijet structure, in Section |V|. We consider the mini- 
jet patterns which arise from additional gluon radiation in the Higgs signal and W + 4 jets 
background, or from 6-quark jets in the ti background, as an additional selection criterion. With 
a central mini-jet veto above prj = 20 GeV, the combined background is reduced well below the 
signal level, without significantly reducing the signal cross section. For an integrated luminosity 
of 100 fb~^, the expected event rate after all cuts is 99 (91) events for a rnn = 800 GeV (1 TeV) 
Higgs boson signal, with a combined background of 41 events. These results suggest that the 
search for the Higgs boson at the CERN LHC can be extended to the 1 TeV region, in the semi- 
leptonic Higgs decay channel. Finally, a summary and our conclusions are given in Section . 



II. CROSS SECTION CALCULATION FOR SIGNAL AND BACKGROUND 

The signal process to be considered at lowest order is the subprocess 

qiq2 ^ qsQAWi^ lu)W{-^ jj) (2) 

and crossing related processes. In the following we require double forward jet tagging (of the jets 
corresponding to quarks qs and ^4) and the presence of at least one additional high transverse 
momentum central jet (from W jj). These requirements are sufficient to eliminate soft and 
collinear divergences and they justify a few approximations in the cross section evaluations which 
will be discussed shortly. 

All cross section calculations are performed numerically, for pp collisions at a center of mass 
energy ^/s = 14 TeV. Individual subprocess cross sections are determined by numerically eval- 
uating polarization amplitudes, mostly by using the amplitude techniques of Ref. |]19[. Even 
though this formalism is well suited to handle massive fermions, all quarks and W^-decay lep- 
tons are treated in the massless approximation, except for the top-quarks. This approximation 
greatly speeds up the calculations. Consistent with it, no Cabibbo-Kobayashi-Maskawa mixing 
is included in the calculation, even for incoming quarks. The error introduced by this approxima- 
tion is well below 5% and, hence, negligible compared to the typical uncertainties of a tree level 
calculation. W decays are evaluated in the zero-width approximation. However, the W — > //' 
decay amplitudes are fully implemented and, thus, all correlations between the decay fermions are 
included in our calculation. Finally, the phase space integrals are performed with the adaptive 
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Monte Carlo integration routine VEGAS pO[ . The statistical error of all Monte Carlo integrals is 
below 1%, except for the W + 4 jets process for which the statistical error on total cross sections 
is ^ 1.5%. 

In all calculations, input parameters are a Z-mass of mz = 91.19 GeV, sin^ 9^ = 0.231 for 
the weak mixing angle, and a = 1/128.75 for the QED fine-structure constant at the electroweak 
scale. From these mw = 79.9 GeV is derived at tree level. The 1-loop formula is used for the 
strong coupling constant as(/i|.), with as(M|) = 0.12. For all processes, the MRS A param- 



eterization of parton distribution functions is used 1 2 1,22]. Even though this parameterization 



is NLO and, hence, we are partially including higher order corrections, these ambiguities in- 
troduce negligible uncertainties. Finally, 6-quark contributions to the initial state are neglected 
throughout. 

In what follows, we give a brief account of calculational details for individual signal and 
background processes. 



A. The Electroweak Process qq — > qq{g)W^W 

The signal process at leading order is qq qqH — > qqW^W~ with subsequent l^-decay, i.e. 
emission of the Higgs boson off a t-channel W or Z a.s shown in Fig. |I](a). For a heavy Higgs boson 
mass {rriH ^ 600 GeV) the narrow Higgs width approximation is no longer applicable and all 
weak boson scattering processes (like the ones shown in Fig. |l|(b)) as well as l^-bremsstrahlung off 
the quark lines (see Fig. |I|(c)) must be considered [|7|, p!3| , p!^ . In principle we need to evaluate the 



full set of 0{aj^) diagrams for a W^W~ final state, including contributions from qq annihilation 
graphs and fermion interchange graphs for identical quarks. We will be requiring a double forward 
jet tag, however, which puts the final state quarks into very different phase space regions and at 
large invariant mass. As a result, annihilation diagrams such as the one shown in Fig. ^(d) and 



the interchange of identical fermions have very small contributions [fq]. They will be neglected 



in the following. Within these approximations the helicity amplitudes for all subprocesses are 



evaluated numerically, using the results of Ref. P 



The signal cross section, as discussed above, contains contributions from non-resonant elec- 
troweak processes such as W bremsstrahlung off the quark lines. Such contributions are in- 
dependent of the mechanism for electroweak symmetry breaking and must be subtracted in 
order to get an estimate of the Higgs boson signal cross section. We model this continuum 
electroweak background by computing the signal at ttlh = 100 GeV. The actual signal cross 
section is then defined as the difference between the heavy Higgs and the ttih = 100 GeV results, 
o'sig = (^{j^h) — (^{mn = 100 GeV). 

In order to understand the characteristics of soft parton emission in the Higgs signal process, 
the 0{aj^as) QCD corrections to the processes of Eq. (H) are needed. The full set of real emission 



diagrams leading to a W^W + 3 parton final state was calculated in Ref. [T^ and we here use 
their results. The subprocesses to be considered are 

qiq2 ^ q^q^g lv)W{^ jj) (3) 

and all crossing related processes like, for example, 

qi9 mmWi^ lu)W{-^ jj) . (4) 
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FIG. 1. Feynman diagrams for the electroweak processes qq — > qqW^W~ . Representative graphs 
are shown for (a) Higgs boson production via weak boson fusion, (b) generic weak boson scattering, (c) 
W bremsstrahlung off the quark hues, and (d) quark-antiquark annihilation. 



Again, s-channel graphs corresponding to gg-annihilation and Pauli interchange graphs for iden- 
tical quarks are neglected. For the Higgs signal calculation at leading order and at 0{as) both 
the renormalization and the factorization scales are set to the smallest transverse momentum of 
the final state partons. 



B. QCD W+ Jets Background 

In signal events with a high transverse momentum W which decays hadronically, W qq, 
the two "jets" in the W decay may merge and form a single high pt jet. In this case the signal 
events produce a W + 3 jets signature. The relevant QCD background for these events comes 
from QCD processes with a W and three jets in the final state. At leading order, two generic 
subprocesses contribute, 



99 qiq29 Iv) 
qiq2 q3q49 W{-^ Iv) 



(5) 



We use the results of Refs. [|T^,^] to calculate the cross sections for these events. All crossing 



related processes are included in the calculation. 

When investigating questions like the ly-invariant mass resolution in ly — > jj decays or the 
additional radiation of soft partons in qq — > qqW^W~ events, the QCD W+ jets backgrounds 
with four partons in the final state are needed. The subprocesses that contribute can be classified 
as 6 quark processes, 4 quark plus 2 gluon processes, and 2 quark plus 4 gluon processes. 
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qiq2 -> qsq^gg lu) (6) 

99 qiq299 lu) . 

The cross sections for these and all crossing related subprocesses were first calculated in Ref. . 
Here we use equivalent matrix elements which were computed by generating the helicity ampli- 
tudes with the program MadGraph IJIB 



For the W^+n jets QCD processes the factorization scale is set equal to the smallest transverse 
momentum of the final state partons. At leading order these cross sections are proportional to 
a", and the strong coupling constant as is evaluated at the corresponding transverse momentum 
of each final state jet, i.e., a" = OILi Q^s(PT,jetJ- 

C. ttjj Background 

For the ti background, the b-quarks from the t —>■ Wb decay are produced mainly in the 
central region, with the two forward jets resulting mainly from QCD radiation. The relevant 
leading order process is the production of ti pairs in association with two jets, which includes 
the following subprocesses 

99 ttgg W^hW~h 99 

qq ttgg W+bW-b gg (7) 
qiq2 ttqiq2 W"^bW~b gig2- 

The exact matrix elements for the O(a^) processes are evaluated, including all the crossing 
related subprocesses. The Pauli interference terms between identical quark flavors in the process 
qiq2 tiqiq2 are neglected, with little effect in the overall cross section rate, due to the large 



differences in the transverse momenta and energies of the final state partons [|T^. The top 
quark decays are simulated in the narrow width approximation, and its mass is set to rrit = 175 
GeV. The structure function scale is chosen to be the smallest transverse energy of the final 
state partons before the top quark decay. The strong coupling constant is evaluated at the 
corresponding transverse energy of the final state partons, prior to the top quark decay, i.e.. 



a 



as{ET{t))as{ET{t))as {pT,iet^ ) as (PTJetj ) • 



In order to study the effects of additional parton radiation in the top quark background, 
one would like to evaluate the tt + 3 jets cross sections as well. Since such a calculation is not 
available yet, we only consider the additional central jet activity arising from the 6-quarks which 
are associated with the top-quark decays. The probability of a 6-quark with > 20 GeV to be 
identified as one of the two forward tagging jets was found to be small (~ 6%) P3]. With the 



transverse momentum and separation requirements on the two tagging jets to be discussed below, 
only this small fraction of the ttjj background is affected by coUinear and infrared singularities. 
Instead of dropping these events altogether we regularize the singularities with the truncated 
shower approximation (TSA). 

D. The Truncated Shower Approximation 

As the transverse momentum of the softest parton becomes small, the perturbative calculation 
of the 0{as) cross section for both signal and background breaks down due to the collinear and 
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infrared divergences associated with gluon emission. In a complete next-to-leading order (NLO) 
calculation these divergences are eliminated due to the cancellation between virtual and real 
emission corrections. For the multi-parton processes considered here, a full NLO treatment is not 
yet possible, however. Instead, we model the effects of multiple soft gluon emission by using the 
truncated shower approximation (TSA) |^. The TSA correctly reproduces the normalization 
of the lowest order cross section (which is free of infrared and coUinear divergences) and it agrees 
with the full NLO calculation when the emission of one additional hard parton is considered. At 
the same time the TSA provides a model for the collective effects of multiple soft parton emission 
in events with n hard jets. Specifically, the tree level n + 1-jet cross section is replaced by 

a[n + 1 ])tSA = K jTSA{PTj,min> -] dpTj,min ■ (8) 

J (lPTj,min 



Here pTj,min is the transverse momentum of the softest jet. 



fTSAiPTj,min) = 1 - exp ( j (9) 

\ Ptsa / 



is a Gaussian cutoff factor, and if is a multiplicative factor that effectively includes the full 1-loop 
corrections. It has been shown that the if-factor for vector boson scattering in pp collisions is 
small {K = 1.06 at the CERN LHC for mn = 800 GeV) |26|. Since fsT-factors are unknown for 
the background processes we set i^' = 1 throughout this study. The parameter ptsa is chosen so 
that the cross section of Eq. (|]) correctly reproduces the lower order n jet cross section. 

For jet transverse momenta pTj < Ptsa, the TSA leads to a reduction of the transverse 
momentum distribution of the hard n-jet system which simulates the canceling of multiple soft 
parton momenta. Thus Ptsa provides an estimate of the jet transverse momentum scale below 
which the emission of multiple soft gluons becomes important. In the phase space regions for 
hard jets to be discussed below, we find values of order Ptsa ~ 40 GeV for the W + 4 jets QCD 
background as compared to Ptsa ~ 8 GeV for the signal. 

For the tt background we apply the TSA only to those events where one of the 6-quarks 
arising from the top quark decays gives rise to at least one of the two forward tagging jets. In 
such events one of the two additional final state partons can be soft, and the cross section is 
enhanced in the region of the phase space where the transverse momentum of that jet is small. 
In order to avoid this singular behavior, the TSA is applied only to this fraction of events, with 



ptt^^ = 42 GeV 1 12 



III. DECIPHERING THE HIGGS SIGNAL FROM BACKGROUND 

The signal process qq — *• qqH —> qq {W — > lu) {W qq) gives rise to two forward tagging jets, 
one (or two) hard central jets from the hadronic W decay, and a leptonic W decay signature. In 
order for the W decay products to be identified, it is required that each event contains a charged 
lepton, / (either e or /x), in the central region, with 

PT/>100GeV, |r7/|<2, ARi, = ^(r^, - r^^f + (0, - 0,)2 > 0.7 , (10) 

where pti is the transverse momentum of the lepton, rji is the lepton pseudo-rapidity, and ARij 
is the distance between the lepton and any identified jet in the pseudo-rapidity-azimuthal angle 
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plane. In addition, it is assumed that each event has large missing transverse momentum due to 
the neutrino of the W ^ lu decay, 

^j, > 100 GeV . (11) 
All final state partons are identified as jets if they satisfy 

PTj> 20 GeV, \rij\<4:.5, (12) 
and if they are well separated in the pseudo-rapidity-azimuthal angle plane, with 



ARjj = V(Ar/)2 + (A0)2 > 0.7 . (13) 

The requirements of Eq. (|12D are superseded by more stringent requirements for the tagging jets 
and for the Higgs decay products. The hadronically decaying W of the Higgs boson signal is 
identified by requiring the existence of a large transverse momentum jet in the central region, 

p^Tj > 300 GeV , \ri^\ < 2 . (14) 

The two quark jets in the process qq — > qqH are tagged by requiring the presence of two 
additional jets, in the forward and backward regions, with 



p^7' > 50 GeV , 2 < r/'"^' < 4.5 , (15) 



tagl 



and 



P?f^>30GeV. {-'i^P^^l (16) 

The asymmetric transverse momentum requirement on the two tagging jets is motivated by the 
fact that one of the jets has substantially higher median p^ than the other, as shown in Fig. 

The resulting cross sections for the signal and the background are shown in the first column 
of Table |I[ The 1^+3 jets background is a factor of 20 larger than the signal, whereas the ti 
background is a factor of three larger. In contrast, the electroweak continuum background is 
strongly suppressed by double tagging, due to the fact that the rapidity distribution of the two 
tagging jets for the electroweak background peaks in the central region [^]. At this level it 
contributes only ~ 10% to the signal cross section. 

For the W+3 jets cross section with the tagging requirements of Eqs. (|T3p and ([T^), it is 
important to ensure that it is a well defined hard scattering process for which a perturbative 
evaluation is reasonable. In order to investigate the effect of double tagging on the cross section. 



we computed the W+1 jet cross section with the single jet satisfying Eq. (p!lD, and the W+2 jets 
cross section with the two jets satisfying Eqs. (|14|) and (^). The W-decaj leptons must satisfy 
the cuts of Eqs. (p!0|JTT|) . The W+1 jet cross section is 2.16 pb, whereas the W+2 jets cross 
section is 0.57 pb. The corresponding reduction factors are 3.8 from W +1 jet to W + 2 jets and 
8.7 from W+2 jets to W+3 jets, respectively. These factors are typical for perturbative QCD 
processes with successively larger numbers of jets and lend credence to the use of perturbation 
theory in the evaluation of the QCD W + n jets backgrounds. 

The signal to background ratio can be further improved by utilizing differences in the tagging 
jet characteristics between signal and background. The two forward jets for the signal are 
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FIG. 2. Transverse momentum distributions (a) for the highest and (b) for the lowest pT forward 
jets at y/s = 14 TeV. For each event, a high transverse momentum lepton in the central region is 
required, pTi > 100 GeV and \r]i\ < 2, as well as missing transverse momentum of > 100 GeV. In 
addition a minimum of three visible jets are required, each with px > 20 GeV. The solid line represents 
the full signal calculation for mn = 800 GeV, while the dotted line represents the continuum electroweak 
background (mjy = 100 GeV). 



very energetic and their energy distributions decline slower than the energy distributions of the 



two forward jets for the W+3 jets and the ttjj backgrounds |l8|,|9|j2^. The softer jet energy 
distributions for the background reflect the fact that these jets tend to come from soft gluon 
radiation in the forward region. By requiring that both tagging jets satisfy 

the signal to background ratio can be improved by more than a factor of two (see second column 
of Table 1). A second distinction arises in the pseudo-rapidity separation of the charged decay 
lepton and the closest tagging jet. For the Higgs signal there is little correlation between the 
two because the leptonic M^-decay arises from the decay of a scalar particle which moves slowly 
in the laboratory frame. By contrast the + 3 jets background contains many events with W 
bremsstrahlung off one of the tagging jets, and such events favor a small separation between the 
jet and the decay lepton. These differences are exploited by imposing a cut, 

- Vi\ > 2 , (18) 

on the separation between the decay lepton and the two tagging jets. The signal and background 
cross sections after all hard cuts are shown in the third column of Table |. The W+3 jets 
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TABLE I. Signal and background cross sections Ba in fb after double jet tagging. The decay lepton 
acceptance requirements are pji > 100 GeV and \rji\ < 2, and each event is required to have missing 
transverse momentum -^j. > 100 GeV. The signal is defined as a{mH) — cr(mH = 100 GeV). 



= 800 GeV 
ruH = 100 GeV 
W+3 jets 
ti + jj 



Double jet 

tagging 
Eq.(|l|)-(|l6|) 
3I5 
0.26 
66.3 
8.01 



signal: 

ruH = 800 GeV 



+jet energy 
cut 
[Eq-dl)] 
1.96 
0.18 
18.2 
3.05 



+ lepton-tagging 
jet separation 
[Eq-dl)] 



1.58 
0.10 
8.36 
1.55 



1.78 



1.48 



background is still a factor of 6 larger than the signal, whereas the ti background has been 
reduced to the same level as the signal. 

We find that any further hardening of the acceptance criteria discussed so far will degrade 
the signal rate appreciably, with only marginal improvement to the signal's statistical signifi- 
cance. Additional information is needed in order to further suppress the background without 
significantly degrading the signal cross section. This is the focus of the following two sections. 



IV. RECONSTRUCTION OF THE W ^ JJ DECAY 

One additional piece of information is provided by the internal structure of the large transverse 
momentum jet in the central region which represents the hadronically decaying W of the Higgs 
signal. The invariant mass of this system, which may or may not be resolvable into two separate 
jets, provides an important criterion for suppressing the QCD H^+jets background. Whenever 
the pair of jets from the hadronic W decay can be resolved, further information is gained. A 
Higgs boson decays mostly into longitudinally polarized Ws whereas backgrounds with real Ws 
are dominated by transversely polarized weak bosons. The angular distribution and the energy 
asymmetry of the two central jets are sensitive to the polarization of the W boson, and can be 
used in order to test whether or not the reconstructed W is the longitudinally polarized decay 
product of the Higgs boson . 



Some of these questions have been studied previously with the aid of parton shower Monte 
Carlo programs p|,p|, |30| , |28[] . Since we have a full QCD matrix element calculation available for 



the production of 1^ + 4 jet events, we can avoid the approximations inherent in a parton shower 
program and use full tree level QCD to simulate the two forward jets as well as the two central 
jets which would fake the hadronically decaying W. The tt background is not included in this 
study of jet mass effects since, similar to the signal process, the observed central jet pair is the 
result of the decay of a real W boson, which, typically, is longitudinally polarized. 

When using W mass reconstruction, the experimental resolution of the dijet mass is the 
limiting factor. In order to model these experimental errors, the lateral granularity of the detector 
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must be taken into account. Following the design specifications of the CMS detector [|3T[, we 
divide the legoplot into cells of size 



Ar/A0 = 0.1 X 0.1 . (19) 

The momentum vectors of the two central jets are then corrected to point to the center of the 
cell. This correction is applied to the Higgs signal and to the continuum electroweak background, 
but not to the W+A jets background. For the former these smearing effects are important, due to 
the resonance in the dijet invariant mass spectrum at rrijj = mw, while the background exhibits 
a fairly fiat dijet mass spectrum which mitigates any smearing corrections. The finite energy 
resolution of electromagnetic and hadronic calorimeters affects both signal and background cross 
sections because energy and transverse momentum distributions are typically quite steep. These 
energy resolution effects are taken into account by Gaussian smearing of the overall energy scales 
of massless parton four- momenta, with relative energy uncertainties [E^ , 



AEem 0.03 0.15 , , 

-^ = -/^®^®0.005, (20) 



and 



AE,a, _ j tI ® ¥ ® 0-03 when |r/,| < 2.5 



E I ^ (g, 3X) ^ ^i^gn 2.5 < IriA < 4.5 

I VE ^ ^ 



(21) 



Here '(8>' means that the terms are added in quadrature. Energy smearing according to 
Eqs. (En,pT|) is applied to both the signal and the W + 4 jets background. 



A. Resolution of Jet Pairs from W — > qq Decay 

The resolution of the two jets from W — > qq decay depends on the angular separation, 

= ^i^n-Vnf + i^n-HY , (22) 

of the two partons in the legoplot. Furthermore, in order to suppress the QCD W + 4 jets 
background, it is advantageous to raise the transverse momentum threshold for each of the two 
central jets above our nominal value of 20 GeV. For the study of W hadronic decay we thus 
require the existence of two central jets, with the following acceptance requirements which are 
added to the requirements of the previous section: 

1. Each of the central jet candidates must have large transverse momentum and be in the 
central rapidity region, 

p'tj > 50 GeV, < 2 . (23) 

All jets passing this criterion form candidate pairs for the hadronic W decay products. 
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2. For each candidate pair, the reconstructed W must have large transverse momentum, and 
it must he in the hemisphere opposite to the lepton-neutrino pair, 

> 300 GeV, -(j)i^\> 90°, (24) 

where (t)jj — ct)iv is the azimuthal angle between the jet-jet pair and the lepton-neutrino pair. 

3. Finally, it is required that the two central jet candidates are separated by 

0.2 < /\m- < 1.0 . (25) 
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FIG. 3. Separation between the two central jets in the pseudo-rapidity-azimuthal angle plane. 
The solid histogram represents the full signal calculation for mn = 800 GeV, which still contains the 
continuum electroweak background (dotted line), calculated in terms of the mn = 100 GeV cross section. 
The QCD W+A jets background is given by the dash-dotted histogram. Energy smearing according to 
Eqs. (2^,21) is applied to signal and background processes. Finite detector granularity (see Eq. (p!s|)) is 
taken into account for the electroweak processes. 



The '^Rjj distributions for the signal and the background are shown in Figure ^. The min- 
imum ^Rjj requirement of 0.2 is still sufficient to eliminate the final state colhnear singularity 
of the W+4 jets cross section at ^Rjj — 0. Notice also that the maximum separation cut, 
^Rjj < 1, has an appreciable effect on the -|- 4 jets cross section only. Due to the large 
transverse momentum carried by the W (see Eq. (|2^)) and the associated strong boost of the W 
decay products, the two quarks from the W decay are rarely separated by more than ^Rjj = 0.7 
in the laboratory frame. 
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TABLE II. Signal and background reduction factors resulting from an analysis of the central cluster 
which is a candidate for hadronic W decay. The first column gives the efficiency of reconstructing two 
jets in the central cluster, within the cuts of Eqs. (p^^-p5|). Requiring the invariant mass of these two 
jets to lie in the Mw i 15 GeV window yields the additional reduction factor of the second column. 
The product of the two yields the total efficiency listed in column three. 





[Eq. ®-(|2|)] 


a(W^n,MnCut) 
[Eq. (H)] 


Efficiency 


mH = 800 GeV 


0.87 


0.83 


0.71 


mH = 100 GeV 


0.67 


0.86 


0.57 


W+jets 


0.27 


0.21 


0.055 


signal: 








= 800 GeV 


0.88 


0.83 


0.73 



In the previous Section, no requirement was imposed on the internal structure of the central 
hard jet. Resolving it into two jets, corresponding to the W — >• qq' decay, will lead to a reduced 
rate for the signal. The corresponding cross section reduction factors, after the cuts of Eqs. (P5|- 
p5|) , are listed in the first column of Table |I[ For the electroweak processes, the effects of detector 
granularity and of energy smearing must be included also in the determination of the single jet 
cross sections which serve to normalize these reduction factors. For the signal cross section, 
86.4% of the events pass the selection criteria. The reduction in the cross section comes mostly 
from the minimum px requirement of 50 GeV for each of the two central jets. Only two thirds 
of the electroweak continuum background survive the cuts. 

The continuum W+ jets background is reduced by approximately a factor 4 when the resolu- 
tion of the central jet into two hard jets of px > 50 GeV is required. This reduction is gratifying 
since it indicates that the use of perturbative QCD is still warranted, in spite of the small minimal 
separation of 0.2 for the two almost collinear partons which mimic the W — jj decay. 

B. Reconstruction of the W Invariant Mass 

A further reduction of the background is achieved by requiring that the two hard central 
jets are consistent in invariant mass with a hadronically decaying VT-boson. The reconstructed 
invariant mass distribution of these two jets is shown in Fig. ^. 

The distributions for the signal and the continuum background are narrow and peak at M^y. 
The distribution for the W+4 jets background, on the other hand, is flat, reflecting the fact 
that in this case the jet pair is not the decay product of a real W boson. These differences are 
exploited by a simple invariant mass cut on the central dijet pair, 

Mw - 15 GeV < M^j < M^y + 15 GeV . (26) 

The chosen mass window of ±15 GeV is motivated by the signal width in Fig. ^ and agrees with 
results of a more complete detector simulation [|T],|2[. The reduction factors for both the signal 
and the background, due to the dijet mass cut of Eq. (p6|), are given in the second column of 
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FIG. 4. Reconstructed dijet invariant mass distribution of the hadronically decaying W. The sohd 
histogram represents the full signal calculation for uih = 800 GeV, with the continuum electroweak 
background {niH = 100 GeV) given by the dotted line. The QCD W+A jets background is given by the 
dash-dotted histogram. Finite detector resolution is taken into account as in Fig. ^. 



Table 0. The W+i jets background is reduced by an additional factor of 5, whereas 83% of 
the signal events survive the cut. The overall efficiency of the central jet resolution and the W 
invariant mass cut is given in the third column of Table For the Higgs boson signal, 73% of 
the events survive all cuts, in contrast to only 5.5% of the events for the W+i jets background. 
At this level, the W+A jets background is a factor of 2.3 smaller than the signal. Even if the 
central jet pair cannot be resolved, it may still be possible to measure the invariant mass of the 
broad central jet representing the hadronically decaying W boson. The reduction factors in the 
second column of Table |I| and the cross section values in the last column of Table | indicate that 
the W invariant mass cut would reduce the W+ jets background to the level of the signal. 



C. Measurement of the W Polarization 



Any polarization of the hadronically decaying W affects the angular distributions of the two 
W decay jets. A transversely polarized W yields a 1 + cos^6'* distribution whereas longitudinally 
polarized Ws produce a sin^6'* distribution. Here 9* is the polar angle of one of the decay jets 
with respect to the W direction, in the W rest frame. The approximate alignment of the thrust 
axis with the W direction for transverse Ws produces two jets of quite different energies after 
boosting into the laboratory frame. Longitudinally polarized Ws, on the other hand lead to 
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FIG. 5. Energy asymmetry distribution of the two central jets, which are identified as the W jj 
decay products. No dijet invariant mass cut is imposed. The sohd histogram represents the fuh signal 
calculation for rriH = 800 GeV, while the continuum electroweak background {rriH = 100 GeV) is given 
by the dotted line. The dash-dotted histogram represents the QCD W+A jets background. Finite 
detector resolution is taken into account as in Fig. 



approximately equal jet energies. The energy asymmetry, A, of the two central jets therefore is 
an excellent variable to confirm the longitudinal polarization of the Ws expected in Higgs boson 
decay 0. It is defined as 

A = L^LZ^ (27) 
El + E2 

where Ei and E2 are the energies of the two central jets in the laboratory frame. The energy 
asymmetry distributions for the ttlh = 800 GeV signal and for the electroweak and the 11^ + 4 jet 
QCD backgrounds are shown in Fig. |, without imposing the dijet invariant mass cut of Eq. (^61) . 
Each distribution is normalized to the corresponding integrated cross section. The difference 
between the distributions for the signal and the background is striking. For the longitudinally 
polarized signal, the two jets have very similar energies so the distribution peaks at A = 0. In 
contrast, the two central jets of the W + 4 jets background and the two jets arising from the 
decay of a transversely polarized W in the electroweak continuum background have substantially 
different energies, with the distribution peaking at large values of A. 

Clearly, an energy asymmetry cut {e.g. at A = 0.5) would further improve the signal to 
background ratio (by a factor of about 1.6). We do not impose such a cut here because there are 
other tools, namely a jet veto on the additional minijet activity in the central region, which can 
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be exploited for an adequate background rejection. The energy asymmetry can then be used to 
confirm the observation of longitudinally polarized Ws from the Higgs boson decay. 



V. CENTRAL JET VETO 

In contrast to the H — ^ W^W~ signal, the two central Ws will be accompanied by close-by 
6-quark jets in the ttj] background and, as we shall see, the W+ jets QCD background also 
produces more observable central jets than the signal process. A veto on any additional central 
jet activity will thus substantially improve the signal to background ratio |[T^JT^ Jl|,3l ■ 



In order to study the effects of semi-soft parton radiation for the Higgs signal and the back- 



ground, we use the TSA of Section [II D|, and thus we first need to estimate the TSA scales, 
Ptsa^ for the various processes. For W + 4 parton production we find that PrsA = 40.5 GeV 
reproduces the + 3 jets cross section of 8.36 pb in Table |. As discussed in Section fl D| we use 
Ptsa ~ 42 GeV for those ttjj events in which a 6-quark arising in a top quark decay produces one 
of the forward tagging jets. For the Higgs signal, a separate estimate for the mn = 800 GeV and 
niH = 100 GeV cases gives two different values for ptsa which, if used, will lead to an incomplete 
subtraction of the continuum electroweak background. Instead, we match the difference between 
the two cross sections, Ba{mH = 800 GeV) — Ba{mH = 100 GeV), to the lower order cross 
section, which gives Ptsa^^'^ = 7.9 GeV. 

The differences in ptsa values between the signal and the background reflect the different 
characteristics of the corresponding hard scattering processes. For the signal, the momentum 
transfer, Q, to the color charges is given by the virtuality of the incident weak bosons in the 
longitudinal weak boson scattering process and, hence, Q signal ~ Pt^ ~ M^. For the W + 4 jets 
and the ti backgrounds, on the other hand, the corresponding scales are substantially larger, of 
the order of Et(W) or, even, the partonic center of mass energy. Very roughly, ptsa, the jet 
transverse momentum scale at which multiple parton emission becomes important, is one to two 
orders of magnitude smaller than the momentum transfer of the corresponding hard scattering 
process. 

The signal and the background cross sections which are obtained within the TSA are given in 
the first column of Table |IT1[ Within the Monte Carlo errors they agree with the 1^-1-3 jets, the 
tijj and the signal cross sections in the last column of Table |. For all results presented in this 
section, the hadronically decaying W is again assumed to decay into a single observable jet. The 
resolution of this jet into two subjets and the effect of an invariant mass cut on this dijet system 
is then taken into account by multiplying with the efficiency factors given in Table |T[ Like 
the Higgs signal, the tt background contains a predominantly longitudinally polarized W which 
decays hadronically and the dijet resolution and dijet mass cut efficiencies for these decays will 
be similar to the ones found for the signal. The signal efficiency of 0.71 has therefore also been 
used for the tt background in the second column of Table |T|. This procedure gives a conservative 
estimate of the top quark background since the combinatorial dilution of the W — > jj peak, due 
to b quarks misidentified as W decay jets, is not taken into account. 

For the Higgs signal the two forward tagging jets define the phase space region in which 
to veto minijet activity. Color coherence favors additional parton emission outside the rapidity 
range bounded by the two tagging jets. A good way to capture the differences between the signal 
and the various backgrounds is by plotting the cross sections as a function of Ar/^j, the smallest 
relative distance, in units of pseudo-rapidity, between the extra jet and the two tagging jets. 
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TABLE III. Signal and background cross sections Ba in fb, before and after the veto of additional 
central jets. Cuts in the first column are the same as in the last column of Table |, but cross sections 
are obtained within the TSA. The second column includes the dijet resolution and mass reconstruction 
efficiencies of Table || and columns three and four give cross sections and expected event rates after the 
central jet veto of Eq. (H). The Higgs signal cross section is defined as Baimn) — Ba{mH = 100 GeV). 







Hard cuts 


~ — r~ 

+ jj resolution 


+ central jet 




Number of events 




+ soft jet 


and Mjj cut 


veto 


£ = 100 fb"^ 




(TSA) 


efficiency 






mH = 800 GeV 


1.64 


1.17 


1.07 


107 


= 100 GeV 


0.17 


0.10 


0.08 


8 


niH = 1 TeV 






0.99 


99 


W+4 jets 


8.49 


0.47 


0.21 


21 


tt 


1.55 


1.10 


0.12 


12 


signal: 










mH = 800 GeV 


1.47 


1.07 


0.99 


99 


ruH = 1 TeV 






0.91 


91 



sign 



?7*"^( closest) 



soft 



(28) 



Here ?7*"^ (closest) is the pseudo-rapidity of the forward tagging jet which is closest to the soft 
jet. The sign in Eq. (|28|) is chosen such that Arjjj is negative if the additional jet is outside 
the pseudo-rapidity interval bounded by the two tagging jets and positive otherwise. The Arjjj 
distribution for additional jets with ptj > 20 GeV is shown in Figure ^. Both background 
distributions peak at Arjjj ~ 2, indicating that the additional jet is predominantly emitted in 
the central region, between the two forward tagging jets. This is in contrast to the signal process, 
where an additional jet is emitted more forward than the tagging jets ~ 50% of the time. Because 
of the small scale governing gluon emission in the signal, only a small fraction of all signal events 
has an extra parton with transverse momentum in excess of 20 GeV. A good strategy is therefore 
to veto events with any additional jets between the two tagging jets, i.e. events which have an 
additional jet satisfying 



soft ^ 

Ptj > PT,veto 



20GeV, 



soft 



tanl taq2 



(29) 



The rapidity requirement corresponds to a cut Arjjj > in Fig. ^ 

The probability for finding a veto jet candidate depends strongly on the minimum transverse 
momentum, pT,veto, of the additional jets. Within the TSA this probability can be estimated by 
integrating daTSA/dPrP ^'^^^ allowed pTj range. The result is then normalized to the lowest 
order cross section. Thus, 



soft 



> Pt, 



veto 



darsA 



O'LO JpT.veto dp' 



soft 
Tj 



1 soft 

dPxj 



(30) 



with r]. 



soft 



tanl taq2 



This probability is shown in Figure |^. At PT,veto = 20 GeV, the proba- 



bility for finding a veto jet candidate in a signal event is below 10%, whereas there is substantial 
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FIG. 6. Rapidity distance, Arjjj, of secondary jets oipTj > 20 GeV from the closest forward tagging 
jet. Results are shown for the mn = 800 GeV 0{as) electroweak processes (solid line), tijj production 
(dash-dotted line) and QCD VF+jets production (dotted line). Negative values of Ar/jj correspond to 
secondary jets outside the rapidity interval formed by the two forward tagging jets. 



probability to find such a jet in background events (~ 55% for the W+A jets and ~ 90% for the 
tijj backgrounds). For the ti background, the veto probabihty tends to 1 as PT,veto 0, due 
to the fact that one of the two 6-jets is almost always emitted in the veto region. Within the 
TSA, the veto probability for the W+A jets background remains substantially less than 1, even 
if PT,veto —>■ 0. This happens because in the TSA only one additional parton is emitted, with 
finite probabihty to be outside the veto region, as seen in Figure ^. Thus, at small PT,veto values, 
the TSA underestimates the veto probability. 

An improved estimate of the veto probability at low PT,veto values is obtained by assuming 
that in the soft region multiple parton emission is dominated by the emission of gluons, and that 
the gluon emission probability exponentiates. This model of multiple minijet emission predicts a 
Poisson distribution for the multiplicity of additional minijets in hard scattering events. Indeed, 
recent CDF data are well described by this ansatz Within this exponentiation model, 

the veto probability can then be estimated as [O] 



exp 



> Pt, 



veto 



exp 



n+l J soft 

7^ ' "-Ptj 



(31) 



where dan+i/dp^p is the unregularized n + l parton cross section, i.e. the higher order cross 
section without the truncated shower approximation. 
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FIG. 7. Probability to find a veto jet candidate above a transverse momentum pT,veto between the 
two forward tagging jets. Results are derived in the TSA for the mn = 800 GeV electroweak signal at 
0{as) (solid line), ttjj production (dotted line) and QCD T^+jets production (dash-dotted line). For 
QCD VF+jets production the result for soft parton exponentiation is shown as the dashed line. See text 
for details. 



A veto probability estimate based on Eq. (0) is also shown in Figure 0, for the W+i jets 
background. At PT,veto = 20 GeV, the TSA underestimates the veto probabihty by 20%. For 
large values of PT,veto {PT,veto > 50 GeV), the two calculations give essentially the same veto 
probability. In the following estimates for the observability of a heavy Higgs signal at the LHC, 
the more conservative TSA results for the veto probability are used. 

Signal and the background cross sections, after applying the central jet veto with pTj > 
PT,veto = 20 GeV, are given in the third column of Table |IT1[ As expected, the jet veto is 
extremely effective in removing the tt background, due to the presence of the two 6-jets. The 
integrated cross section is reduced by one order of magnitude and is now only ~ 12% of the 
signal cross section. The W+A jets background is reduced by a factor 2.2. These background 
reductions are achieved with a very high efficiency for retaining the signal, with approximately 
91% of the signal events passing the veto criterion. 

After the veto of additional jets in the central region, the signal cross section rate is a 
factor of 2.5 larger than the combined background rate. Assuming an integrated luminosity 
of £ = 100 fb~^, the expected number of events for the Higgs signal and for the background are 
given in the last column of Table |IT1[ At tuh = 800 GeV, 99 signal events are expected with a 
total background of 41 events. These numbers indicate that the Higgs boson can be discovered 



19 



in the H WW lujj decay mode, for Higgs masses up to tuh = 1 TeV. 



VI. DISCUSSION AND CONCLUSIONS 

Our analysis of the H WW — lujj decay mode of a heavy Higgs boson is based on 
complete tree level QCD calculations of the cross sections for the qq qqWW signal as well 
as for the W + 3,4 jets and tijj background processes. Full QCD matrix elements provide 
the most reliable predictions for event features like hard jet distributions, the momentum scales 
governing the emission probability of additional soft jets, and the angular distributions of such 
additional jets. With QCD matrix elements for {W —>■ li'){W —>■ jj) + 2,3 jet events for the 
signal, + 3, 4 jet events for the QCD background and {t — > Wb) (i — > Wb) + 2 jet production for 
the top-quark background we have analyzed optimal criteria for double forward jet tagging, the 
expected resolution and expected background suppression when searching for aW —>■ jj invariant 
mass peak, prospects for measuring the longitudinal polarization of the hadronically decaying 
W of the Higgs boson signal, and we have studied momentum scales and angular distributions 
of additional soft jet emission which would be affected by a central jet veto. 

Previous analyses by the ATLAS 0] and CMS [|l| Collaborations used parton shower programs 
like PYTHIA [03 instead, which give a more detailed description of other aspects of signal 



and background events, like particle content, higher soft jet multiplicities, and the presence 
of an underlying event. Carrying the simulation to the particle level also allows for a more 
realistic assessment of detector response. With fairly similar acceptance cuts on the Higgs decay 
products and using double forward jet tagging and central jet vetoing techniques these studies 
arrived at qualitatively the same answer: that a heavy Higgs boson can be discovered in the 
H WW Ivjj decay mode. It is reassuring that also quantitatively the agreement is 
excellent. The predicted rates for the Higgs boson signal and the various backgrounds found in 
the ATLAS and CMS analyses are somewhat smaller than ours, by up to a factor two. While 
this general trend is expected from including detector efficiencies, our central-jet- vetoing and 
forward-jet-tagging criteria differ sufficiently from those used in Refs. that this agreement 
may be fortuitous to some extent. Since we do not have the parton level cross sections for these 
analyses, a more explicit and quantitative comparison with the ATLAS and CMS studies is not 
feasible. A few differences are noteworthy, however. 

Both the CMS and the ATLAS analysis find suppression of the W+ jets background by a 
factor ^ 3 from a central jet veto, but with very different values of PT,veto = 40 GeV (CMS) and 
PT,veto = 15 GeV (ATLAS), and ATLAS reports a decrease to a factor 2.5 when pT,veto = 40 GeV 
is used, at high luminosity. Our analysis indicates that the veto probability should vary more 
strongly with pT,veto (see Fig. |^) and this question deserves further study. Another difference is the 
transverse energy requirement for the two tagging jets: our choice of transverse energy threshold 
for the tagging jets, ETjitag) > 50 (30) GeV, is motivated by the signal distributions in Fig. |[ 
It is considerably harder than the values Erjitag) > 15 GeV (ATLAS) and Exjitag) > 10 GeV 
(CMS) used in the Technical Proposals of the two detectors [|l],@. Thus our choice definitely is 
conservative. 

Neither of the presently available analyses will be definitive by the time the LHC experiments 
start taking data. The tools presented here, based on state-of-the-art QCD matrix elements, 
can be used to calibrate parton shower Monte Carlo programs. Combined with full detector 
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simulations they will provide the more reliable predictions for longitudinal weak boson scattering 
signals and background processes which are needed to understand LHC data. 
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